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We examine a nonlinear magnetoinductive dimer and compute its linear and nonlinear symmetric, 
antisymmetric and asymmetric modes in closed-form, in the rotating-wave approximation. A linear 
stability analysis of these modes reveals that the asymmetric mode is always stable, for any allowed 
value of the coupling parameter and for both, hard and soft nonlinearity. A numerical computation 
of the dimer dynamics reveals a magnetic energy selftrapping whose threshold increases for increasing 
dimer coupling. 



Metamaterials are novel artificial materials character- 
ized for having unusual electromagnetic wave propaga- 
tion properties, such as a negative dielectric permittiv- 
ity and negative magnetic permeability over a finite fre- 
quency range. This feature makes them attractive for use 
as a constituent in negative refraction index materials 
A subclass of those metamaterials, the magnetic meta- 
materials (MMs), exhibit significant magnetic properties 
and negative magnetic response up to terahertz and op- 
tical frequencies [2|, [3[. 

One of the most well-known MMs consist of a metal- 
lic composite structure consisting of arrays of split-ring 
resonators (SRRs). The theoretical treatment of such 
structures relies mainly on the effective-medium approx- 
imation where the composite is treated as a homogeneous 
and isotropic medium, characterized by effective macro- 
scopic parameters. The approach is valid, as long as the 
wavelength of the electromagnetic field is much larger 
than the linear dimensions of the MM constituents (j-[6| . 

The shortest array is the dimer, and consists on two 
SRRs coupled inductively. In spite of its simplicity, 
a dimer is capable of rich phenomenology including 
magnetic energy transfer, mode localization and even 
chaos Q. Linear magnetic dimers have been used as 
constituent units for envisioned three-dimensional meta- 
materials or 'stereometamaterials' [8[. The dynamics of 
magnetic energy localization in an asymmetric nonlinear 
dimer with dissipation and driving has been examined 
in ref.[7], where it was concluded that asymmetry gives 
rise to chaotic dynamics and is also necessary for strong 
localization in one of the dimer sites. Also, the linear 
magnetic dimer with gain/loss terms constitutes yet an- 
other finite system with PT-symmetry properties, pos- 
sessing a parameter window inside which its dynamics is 
bounded Q. 

In this work we focus on a nonlinear dimer system, 
consisting of two identical SRRs, and examine the ex- 
change and localization of magnetic energy via a closed- 
form computation of the linear and nonlinear modes. We 
also examine the dynamics of selftrapping of magnetic 
energy as a function of the mutual coupling between the 
rings. 
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Figure 1: Magnetoinductive dimer. In (a) M < 0, while in 
(b) M > 0. 



Let us consider an array consisting of two identical split 
ring resonators (SRRs), as shown in FigJU The system is 
characterized by the values of the self-inductance L, lin- 
ear capacitance C/, mutual inductance M, characteristic 
voltage U c and linear dielectric constant e\. In the pres- 
ence of dissipation and driving, the coupled equations 
for the time evolution of the charges Qi,Q 2 are given 
approximately by 
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The most commonly used dimer configurations are in-a- 
plane and on-an-axis (FigJTJ. In the first case, the mu- 
tual inductance is negative, while on the second case it 
is positive. 

Linear case. We start by examining the simplest case, 
where there is no dissipation and driving force, and 
nonlinear effects are negligible. With the definitions 
Q c CiU c ,ujQ \/>/LCi,t = uj t,\ = M/L < 1 and 
q = Q/Qc, Eqs.(pQ) become: 
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Figure 2: Evolution of SRR energies for two coupling values: 
8/10 (left) and 9/10 (right). The solid (dashed) curve denotes 
Hi(H2). In the first case, the evolution is periodic; in the 
second it is quasiperiodic. 



Eqs.(|2|) become 
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where x = ol/%€i. In order to have an oscillatory behav- 
ior at all, we need that q 2 2 < 1/x f° r X > 0- For \ < 
no such restriction is necessary, since the potential is al- 
ways hard in that case. We see that the change A — >• — A 
is equivalent to qi -)> gi,g 2 -> ~Q2 or q 1 -> -<?i,<?2 -> q2- 
We could call this the "staggered" mode. More interest- 
ingly, Eq.([7j) admits the conserved quantity: 



It is easy to see that H = (1/2) (g^ + gfHt 1 / 2 ) (<?i+<?2) + 
\q1q2 is a constant of motion, dH/dt = 0. We define the 
energy contents Hi and H2 of each SRR as 
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We now look for stationary solutions: qi^ ~ 
exp(i/3r). After inserting into (|2j), one obtains j3 = 
±1/Vl + A, ±1/Vl — A. As a direct application, we 
consider the initial- value problem <?i(0) = <7o?(h(0) = 
0?^2(0) = 0,^2(0) = 0, whose solution is 
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Since the two frequencies involved, 1 / yl ± A are in gen- 
eral incommensurable, the motion of qi{r) and ^(t) 
will be quasiperiodic, unless A can be written as (p 2 — 
q 2 )/(p 2 + <? 2 ), where g are integers. Once in posses- 
sion of q\(r) and ^(t) in closed form, we could write 
Hi and H2 explicitly as function of time, but the expres- 
sions are rather cumbersome and not particularly illu- 
minating. Figure [2] shows examples of the evolution of 
the SRRs energies for similar coupling values that lead 
to qualitatively different evolutions. Besides the typical 
presence of quasiperiodicity, the SRRs energies are not 
positive-definite. We see that even in this simple case, 
the dynamics is considerably more complex than for the 
linear case of the Discrete Nonlinear Schrodinger (DNLS) 
dimer. In particular, it is not clear now how to define a 
proper coupling time. This is due to the fact that en- 
ergy is not only stored as charges in the capacitors of the 
SRRs, but also in the magnetic fields across the SRRs 
slits. 

Nonlinear case. In the presence of nonlinear effects, 



X, 
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which we call the magnetic energy. The individual energy 
contents of each ring are then 
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Rotating-wave approximation (RWA). We look for sta- 
tionary modes of the form qi(r) = qi sin(Qr), ^(t) = 
q2 sin(^r), and use the approximation sin(x) 3 « 
(3/4) sin(x). We obtain the coupled equations 
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which is invariant under the change qi — >• ^2 and vicev- 
ersa. Also if (#1,^2) is a solution, so will (— <2i, — (72)- Let 
us examine some of the modes implied by Eqs. (p"0|) : 
(i) qi — q2 = q (Symmetric mode). This leads to 
(1 - ^ 2 (1 + X))q - (3/4)x q S = which implies 
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(ii) q2 = —qi = q (Antisymmetric mode). This leads to 
(1 - ft 2 (l - X))q - (3/4)x q 3 = 0, which implies 
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and n 2 < 1/(1 - A) if x > 0; otherwise ft 2 > 1/(1 - A) if 
X<0. 

(iii) q 2 7^ q 2 (Asymmetric mode). After multiplying the 
first eq. in ([TP]) by qi and the second eq. by q2 , and after 
subtracting, we obtain 
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Figure 3: Nonlinear mode amplitude as a function of mode 
frequency, for the symmetric, antisymmetric and asymmet- 
ric modes. Solid (dashed) portions denote stable (unstable) 
regimes. The solid curves labelled gi, q2 refer to the asym- 
metric mode, (x = 1, A = 0.33) 



Next, we multiply the first eq. in (fT0|) by q 2 and the 
second one by q±, we obtain (q 2 — q 2 )[\Q 2 — (3/4)x<Zi#2] = 
0. Then, assuming q 2 ^ q 2 , we obtain 



From Eqs. (fT3|) and (fl4]h we finally obtain 



91 =2 



<?2= 2 



n 2 + 2A^ 2 

"(374k 



ft 2 + 2Aft 2 

(3/4)x 




Q 2 - 2\9? 
(3/4)x 



1 - ft 2 - 2Aft 2 



(14) 



(15) 



(16) 



where, without loss of generality, we have assumed q\ > 
q2. In order to have well-defined real solutions, we require 
|A| < 1/2 and 

Q < Min{ 
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n > Max{ 
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Stability of the RWA modes. Let us examine the lin- 
ear stability of the RWA modes we just found. We set 
Qi( T ) = gisin(ftr) + Sq 1 (r) and q 2 (r) = q 2 sin(^r) + 
(5^2 (t), where |&/i,2| <^ #1,2 • After inserting this into 
Eq.([7j), keeping only linear terms in Sqi^, and make use 
of the RWA for sinx 2 ~ 1/2, we obtain the linear system: 
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Figure 4: Time evolution of individual SRR energies for 
increasing values of initial energies on one of the SRRs. 
if i(0) = 0.73 (left), ifi(0) = 1.0 (right), #i(0) = 1.2 (bot- 
tom). The solid (dashed) curve denotes H\(H2). Inset shows 
small amplitude oscillations (% = 1/6, A = 0.1) 
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We pose Sqi(r) = SAsm(ur) and Sq 2 (r) = 5Bs'm(ujr). 
The system has a nontrivial solution provided 
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where a\^ 2 = 1 — (3/2)x<? 2 2 - From (f2T|) we obtain 
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Now, in order for the modes to be stable, we need uj 2 > 0. 
Case x < 0. In this case, a± = 1 + (2/3|x|) > and 
a 2 = 1 + (2/3|x|) > 0. The condition uj 2 > 0, leads to 
(a\ — a 2 ) 2 > — 4A 2 ceia2, which is always satisfied. Thus, 
for x < all modes are stable. 
Case x > 0- I n this case, we need 



Qi < Q2 < T~ and 
3x 3x 
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Let us apply these conditions to the modes we found 
before. 

(i) qi = q 2 = y/(l-n*(l + \)/((3/4)x). Condition 



(|23j) implies Q 2 > 1/2(1 + A), for stability. Combin- 
ing this stability condition with the existence condi- 
tions, we conclude that the mode exists and is stable 

if 1/(2(1 + A)) < Q 2 < 1/(1 + A). 

(ii) q2 = -qi = Vi 1 ~ « 2 (1 ~ A)/((3/4) X ). Condition 
(|23|) implies Vt 2 > 1/2(1 - A). Combining this stabil- 
ity condition with the existence conditions, we conclude 




2 4 6 

INITIAL ENERGY 




curvatures: Positive for hard nonlinear ity, while negative 
for soft nonlinearity. 
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Figure 5: Long-time averaged fraction of initial energy re- 
siding on SRR as a function of its initial energy, for several 
coupling parameter values : A = 0.05 (solid), A = 0.1 (dotted) 
and A = 0.15 (dashed). Left: Hard nonlinearity (% = —1/6). 
Right: Soft nonlinearity (% = 1/6). 



that the mode exists and is stable if 1/(2(1 — A)) < Q 2 < 
1/(1 -A). 

(hi) q\ ^ q%. Using (fT5j) . (fl6j) and (|23]), one obtains the 
condition Q 2 > 0. Therefore, this mode is always stable 
in its existence domain, at least in the RWA framework. 

Dynamics. An interesting question concerning the dy- 
namics of the dimer is the selftraping problem: If we start 
with an initial condition where all the energy is residing 
on a single SRR, what are the conditions needed for that 
energy to remain "localized" or "selftrapped" on the ini- 
tial SRR. An initial condition of the type q±(0) = qo, 
while (72 (0) = gi(0) = g 2 (0) = 0, ensures iJi(0) = H, 
#2(0) = 0. Figure [H shows the time evolution (from 
Eq.0) of the SRRs energies for three different values 
of initial energies, revealing a clear selftrapping transi- 
tion, for initial energies above some threshold. Now, even 
though the curves resemble the ones observed for the cu- 
bic DNLS dimer, a close examination reveals the presence 
of small amplitude oscillations (not readily apparent in 
FigHJ) around the main tendency. The initial- value prob- 
lem seems, in fact, more difficult than its DNLS counter- 
part, and the possibility of achieving a closed- form solu- 
tion seems doubtful at this stage, unless we find a way 
to average over the extra frequencies. 

Figure [5] shows the long-time averaged fraction of en- 
ergy remaining on the initial SRR, as a function of its 
initial energy, for hard and soft nonlinear it ies, and for dif- 
ferent mutual inductance couplings. Clearly, for a given 
value of the nonlinearity parameter, the threshold en- 
ergy needed for selftrapping is an increasing function of 
the inductive coupling A. This information is conveyed in 
a more clear manner in Fig |6] where we show the thresh- 
old value of Hi (0) = H needed to effect selftrapping as a 
function of the inductive coupling parameter. For both 
cases, the "soft" and the "hard" nonlinearity cases, the 
minimum initial energy is a increasing function of cou- 
pling, as expected on general grounds, but with different 
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Figure 6: Minimum initial energy on a SRR to begin selftrap- 
ing, as a function of the value of inductive coupling. Left: 
Hard nonlinearity (% = —1/6). Right: Soft nonlinearity 
( X = l/6). 

Conclusions. We have computed in closed form the 
stationary modes of a nonlinear magnetoinductive dimer 
in the rotating- wave approximation. The linear stability 
window for the symmetric, antisymmetric and asymmet- 
ric modes was obtained in closed form, finding that the 
asymmetric mode is always stable. The dynamics evolu- 
tion of a localized initial excitation reveals a selftrapping 
transition, with an energy threshold that increases with 
an increase in magnetic coupling. 
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